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ABSTRACT 
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We present the results of six-epoch Very Long Baseline Array (VLBA) observations of 3C 66A. The high-resolution Very Long 
Baseline Interferometer (VLBI) maps obtained at multi-frequency (2.3, 8.4, and 22.2 GHz) simultaneously enabled us to identify 
the brightest compact component with the core. We find that the spectrum of the core can be reasonably fitted by the synchrotron 
self-absorption model. Our VLBA maps show that the jet of 3C 66A has two bendings at about 1.2 and 4 mas from the core. We also 
give possible identifications of our jet components with the components in previous VLBA observations by analysing their proper 
motions. We find consistent differences of the position from the core in one component between different frequencies at six epochs. 

Key words, galaxies: jets - quasars: individual (3C 66A) - radio continuum: galaxies 
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1. Introduction 

The source 3C 66A (B0219+428, 4C 42.07) is a low-frequency 
peaked BL Lac object (LBL) peaking at IR-UV wavebands (e.g., 
Petri et al. |2003t Joshi & Bottcher |2TjQ6l ). Since its optical coun- 
terpart was identified by Wills & Wills ( I1974| l, a large number 
of observations from radio to gamma-ray bands have been per- 
formed on 3C 66A. As a blazar, 3C 66A also exhibits a promi- 
nent variability at radio, IR, and optical (cf. Bottcher et al. |2005l 
hereafter B05). 

Price et al. ( I1993I ) gave the first Very Long Array (VLA) 
maps at 1.5 and 5 GHz, both showing that 3C 66A has a 6 
extended structure along the position angle (PA.) ~ 170°. New 
VLA maps at 1 .5 and 5 GHz (Taylor et al. |1996l l further revealed 
that 3C 66 A has two weak lobes, one at 15 south of the central 
region, the other at 10" with PA. about -20°. Taylor et al. ( 119961 
present the first 5 GHz VLBA map, exhibiting a typical core -jet 
structure in 3C 66A. Jorstad et al. (120011 hereafter J01) stud- 
ied the jet's kinematics of 3C 66A through multi-epoch VLBA 
observations at 22 and 43 GHz, and they detected superlumi- 
nal motions for four jet components with the apparent velocities 
from 15.3c to 29.7c (the apparent velocities shown in this pa- 
per have been calculated by using the cosmological parameters 
given below), and one stationary jet component. A more detailed 
kinematic study of 3C 66A was done by Jorstad et al. ( 120051 
hereafter J05) through the total and polarized intensity VLBA 
observations at 43 GHz at 17 epochs spanning 1998 March to 
2001 April, with some epochs accompanied by nearly simulta- 
neous polarization measurements at 100, 222/353 GHz, and the 
optical wavebands. B05 presented a coordinated observation of 



3C 66A during the Whole Earth Blazar Telescope (WEBT) cam- 
paign at multiple frequencies including radio, infrared, optical, 
X-rays, and gamma-rays at very high energy. Their VLBA obser- 
vations at three epochs at 22 and 43 GHz revealed superluminal 
motion for only one jet component with an apparent velocity of 
12.1c, other jet components remain stationary, and the radial ra- 
dio brightness profile suggests a magnetic field decay oc r _1 and, 
thus, a predominantly perpendicular magnetic field orientation. 



We describe in Sect. 2 our VLBA observations of 3C 66A at 
relatively low frequencies of 2.3 and 8.4 GHz during six epochs 
and at 22.2 GHz during the first two of six epochs. The results are 
presented in Sect. 3. We analyse the spectra of the VLBI com- 
ponents in Sect. 3.1, study the kinematics of the jet components 
in Sect. 3.2, and the positional difference of the jet components 
caused by the frequency-dependent core position offset in Sect. 
3.3. The discussion of results is given in Sect. 4, followed by 
conclusions in Sect. 5. 



Send offprint requests to: Z.-Q. Shen, e-mail: zshen@shao .ac.cn 



We adopt 0.444 as the redshift of 3C 66A (Miller et al. |T978j >, 
although Bramel et al. (2005) thought that this redshift value 
is still uncertain and deserves more spectroscopic observations. 
Throughout this paper the radio spectral index a is defined as 
S v oc V, By assuming Hq = 71 km s _1 Mpc~', Q.m - 0.27, and 
Q. A = 0.73 (Spergel et al. 12003]) , we have a scale for 3C 66A (z 
= 0.444) of 1 mas = 5.69 pc. 
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2. Observations and data reduction 

The observations of 3C 66A were performed with the NRACQ 
VLBA at 2.3 and 8.4 GHz during the following six epochs: 
2001 March 13-14 (2001.20), 2001 June 25 (2001.48), 2001 
November 9 (2001.86), 2002 February 8-9 (2002.11), 2002 
February 21-22 (2002.14), and 2002 June 14 (2002.45). The 
VLBA observations at 22.2 GHz were successfully performed 
at the above first two epochs. These observations were used to 
search for the binary black hole in 3C 66B (Sudou et al. I2003t . 
The on-source time at each epoch is about 30 min at 2.3 and 
8.4 GHz, and 70 min at 22.2 GHz. The recording bandwidth is 
16 MHz (two 8 MHz IF channels) at 2.3 and 8.4 GHz, and 32 
MHz (four 8 MHz IF channels) at 22.2 GHz. All the data were 
recorded with the right-circular polarization (RCP) mode at 2.3 
and 8.4 GHz and with the left-circular polarization (LCP) mode 
at 22.2 GHz. 

The data correlation was done with the VLBA correlator in 
Socorro, New Mexico, USA. Post-correlation data analysis was 
done in the NRAO AIPS software (Schwab & Cotton [TW3l for 
fringe fitting and in the Caltech DIFMAP package (Shepherd 
1 1 997b for hybrid mapping. A prior visibility amplitude calibra- 
tion was done using the antenna gain and the system tempera- 
ture measured at each station. During the process of amplitude 
calibration, the correction for atmospheric opacity was also ap- 
plied at all the three frequencies. During the hybrid mapping, 
the CLEAN and phase-only self-calibration were iteratively used 
for the early processing, and the amplitude self-calibration was 
added for the later processing. The final VLBA maps are shown 
in Fig.Q] The mapping parameters of Fig. Q] are listed in Table Q] 
The quantitative description of the source structure was deter- 
mined by the circular Gaussian model fitting to the calibrated 
visibility data, and the results are listed in Table [2] 

For the VLBI spectral analysis, an accurate absolute calibra- 
tion of the flux density is necessary. Unfortunately due to the 
lack of the corresponding VLBA observations for other com- 
pact amplitude calibrators, we could not perform any absolute 
amplitude calibration. Based on our experience of the study on 
PKS 0528+134 and 3C 138 by the mu lti-freq uency VLBA ob- 
servations (Cai et al. 120061 Shen et al. 120051 ). the errors in the 
absolute flux-density calibration are about 10% at 2.3 and 8.4 
GHz, and about 20% at 22.2 GHz. We also used the Difwrap 
program (Lovell 12000b to estimate the errors of the Gaussian 
model-fitting parameters. The errors of the flux densities in the 
Gaussian model fitting estimated by Difwrap are much smaller 
than the errors in the absolute flux-density calibration, so we 
only adopted the errors in the absolute flux-density calibration 
as the error bars of flux density, and adopted the errors of other 
model-fitting parameters estimated by Difwrap. All the errors 
are listed in Table [2] 



3. Results 

The fitted VLBI components listed in Table [2] are also labelled 
in Fig. Q] We label the component at the north end (having the 
smallest size and the highest brightness temperature) k compo- 
nent, which is further identified with the core according to its 
spectral characteristics in Sect. 3.1. Other jet components are la- 
belled alphabetically according to their separations from the k 
component. 



1 The National Radio Astronomy Observatory (NRAO) is operated 
by Associated Universities Inc., under cooperative agreement with the 
National Science Foundation. 



Table 3. The spectral indexes of the VLBI components. 









component 






epoch 


k 




d 


e 


f 


g 


2001.20 


0.067 


-1.28+^ 


-0.44±0.09 


-0.52 


-0.11 


-0.54 


2001.48 


0.22 


-4.52±(> 0.32) 


-0.47+0.09 


-0.28 





-2.27±(> 0.47) 


2001.86 


0.052 




-0.51 


-0.11 


-0.48 


-0.19 


2002.11 


0.05 




-0.57 


-0.24 


-0.23 


-0.54 


2002.14 


0.21 




-0.56 


-0.15 


-0.32 


-0.11 


2002.45 


0.03 




-0.41 


-0.12 


-0.17 


-0.30 



Note: The errors, if not shown in the table, are all 0.16. 



Component c detected at 22.2 GHz was not detected at 2.3 
and 8.4 GHz, which may be due to the lower resolution at 2.3 
and 8.4 GHz and its relative weakness at 2.3 and 8.4 GHz when 
compared with components b and d. The maps at 8.4 and 22.2 
GHz show that there are two bendings in the jet emission, at 
about 1.2 and 4 mas from component k, with an extended jet 
structure toward PA. ~ 170° at 20-30 mas. This extension to- 
ward the south is consistent with that of the VLA maps at 1.5 
and 5 GHz (Price et al. [19931 T aylor e t al. [T996j > and the VLBA 
map at 5 GHz in Taylor et al. d 19961 1. But the structure toward 
the north that has been detected in the VLA maps at 1.5 and 5 
GHz is not shown in our VLBA maps and the VLBA map at 5 
GHz either (Taylor et al. |19961 >. This may be due to the lobe in 
the north (PA. ~ -20°), detected by VLA at 1 .5 and 5 GHz, being 
too weak and dispersive, and thus resolved by the VLBA. In the 
following subsections, we study the spectra, motions, and the 
positional offsets of these jet components based on our VLBA 
observations. 

3.1. Spectra of VLBI components and identification of the 
core 

The detection of components (k, b, d, e, f, and g) at more than 
one frequency allows spectral analysis of these VLBI compo- 
nents. We plot the spectrum of k component in Fig. [2] (a) and 
(b), and the other five components' spectra in Fig.|2](c). The data 
in Fig. [2] (a), (b), and (c) are from our observations at 2001.20, 
2001 .48, and 2001 .20, respectively. The spectral indexes of com- 
ponent d at the first two epochs (2001.20 and 2001.48) were ob- 
tained from the linear regression on the data points at 2.3, 8.4, 
and 22.2 GHz. The spectral indexes of the other components are 
simply from the two-point measurements at 2.3 and 8.4 GHz, 
except for component b at the first two epochs whose spectral in- 
dexes come from two-point measurements at 8.4 and 22.2 GHz. 
The results from these spectral fits during six epochs are given 
in Table [3j where we can see that the k component's spectral in- 
dexes, different from other components, are positive, implying 
the existence of an absorption. Its inverse spectrum at the low 
frequency in Fig. [2] (b) also confirms this. Combined with its 
compactness, we identify k as the core of 3C 66A. 

Taylor et al. ( 19961 report on the VLBI observations of 
3C 66A at 5 GHz at 1995.65 with a core flux density 0.67 Jy. 
According to the single-dish monitoring by the University of 
Michigan Radio Observatory (UMRAO), the total flux densities 
of 3C 66A at 5 GHz at epochs 2001.20 and 2001.48 of 2.25 
and 2.37 Jy, respectively, are almost the same as that at epoch 
1995.65 of 2.38 Jy (Taylor et al. [Y996I By assuming that the 
same total flux density infers the same core activity, we can es- 
timate the flux density of VLBI core to be 0.67 Jy at 5 GHz at 
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Fig. 1. The naturally-weighted VLBA maps of 3C 66A at 2.3, 8.4, and 22.2 GHz from top to bottom. The circles superimposed on 
the maps represent the Gaussian model components listed in Table 2, whose names are labelled. The crosses represent the restoring 
beams in Table 1 . 
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Table 1. Description of VLBA maps of 3C 66A shown in Fig. Q] 



Restoring Beam 



V 


Epoch 


S peak 


Major 


Minor 


P.A. 


Contours 


(GHz) 


(yr) 


(Jy/beam) 


(mas) 


(mas) 


(deg) 


(mjy/beam) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


2.3 GHz 


9001 70 


fi78 
u.u / o 


J . 1 o 


1 67 


-1.21 






2001 4R 


0.641 


5.96 


3.68 


-1.24 


1 88x("-1 1 2 4 8 1 6 3? 64 1 28 756^ 




2001.86 


0.748 


5.39 


3.87 


-5.16 


1 .70x(- 1 , 1 ,2,4,8, 16,32,64, 128,256) 




200? 1 1 


0.718 


5.4 


3.9 


-7.06 


1 4?v(-l 1 ? 4 8 1 6 32 64 1 28 256^ 




2002.14 


0.598 


5.47 


3.87 


-7.1 


1 .79x(- 1 , 1 ,2,4,8, 16,32,64, 128,256) 




2002.45 


0.584 


5.33 


3.98 


-3.55 


1 .72x(- 1 , 1 ,2,4,8, 16,32,64, 128,256) 


8.4 GHz 


2001.20 


0.653 


1.6 


1.04 


-0.34 


1 .74x(- 1 , 1 ,2,4,8, 16,32,64, 128,256) 




2001.48 


0.719 


1.61 


1.01 


-2.06 


2.34x(-l, 1,2,4,8,16,32,64,128,256) 




2001.86 


0.696 


1.44 


1.08 


-4.23 


2.18x(-l, 1,2,4,8,16,32,64,128,256) 




2002.11 


0.623 


1.46 


1.09 


-7.81 


2.27x(-l, 1,2,4,8,16,32,64,128,256) 




2002.14 


0.586 


1.43 


1.07 


-6.05 


2.28x(-l, 1,2,4,8,16,32,64,128,256) 




2002.45 


0.494 


1.42 


1.06 


-4.22 


2.63x(-l, 1,2,4,8,16,32,64,128) 


22.2 GHz 


2001.20 


0.718 


0.64 


0.41 


0.57 


2.73x(-l, 1,2,4,8,16,32,64,128,256) 




2001.48 


0.638 


0.71 


0.38 


-0.80 


3.50x(-l, 1,2,4,8,16,32,64,128) 



Notes:(l) Observing frequency; (2) Observing epoch; (3) Peak flux density; (4), (5), (6) Parameters of the 
restoring Gaussian beam: the full width at half maximum (FWHM) of the major and minor axes and the position 
angle (P.A.) of the major axis. (7) Contour levels of the map. The lowest contour level is three times the rms 
noise in the maps. 
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Fig. 2. The spectra of VLBI components in 3C 66A: (a) and (b) are the spectrum of core (k component) at 2001.20 and 2001.48, 
and (c) is the spectra of the jet components b, d, e, f, and g at 2001.20. Shown in (d) is an SSA fitting curve to the spectrum of k 
component with all the data points from our observations at 2001.48 except 5 GHz data from Taylor et al. (1996t. 
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2001.48 (as plotted in Fig. 2 (d)). Then we use the synchrotron- 
self absorption (SSA) formulary = Sqv 25 [1 -exp(-r J v°' _2 - 5 )] to 
fit the data at 2.3, 5, 8.4, and 22.2 GHz, where S o is the flux den- 
sity in Jy at 1 GHz when the SSA optical depth at 1 GHz t s » 1 
and a is the optically thin spectral index. The best-fit a = -0.08, 
So = 0.10 Jy, and t s = 7.92, the fitting curve is plotted in Fig. 2 
(d). From the fitting curve we estimate that the peak flux density 
S m is ~0.66 Jy at the inverse frequency v,„ = 6.52 GHz. Using 
the relation ^ ~ 8(^)i(fe)j(^-)-j(i + z )i for the SSA 
model (Kellermann & Pauliny [l981| ), we estimate the magnetic 
field in the core region ~ 1 . 1 mG, here 9, the angular diameter of 
the core of 0.21 mas is the geometrical mean of the angular di- 
ameters at 2.3 and 8.4 GHz. For comparison, the angular size of 
3C 66A at 5 GHz is 0.2 mas from previous VLBA observations 
(Taylor et al. [19961 



3.2. Proper motion and registration of jet components 

J01, J05, and B05 have studied the jet kinematics in detail, and 
they showed that the kinematics of 3C 66A is complicated; some 
components show superluminal motions, some show the inward 
(apparently toward the core) motions, and some show the zero 
proper motions. We also linearly fit the separations of the jet 
components at 2.3 and 8.4 GHz (Fig.0, where the fitting results 
are listed in Table [4] Since the observations at 22.2 GHz were 
at two short-separated epochs, we do not fit the separations of 
the jet components at 22.2 GHz. We can see from Fig. [3] and 
Table|4]that there are no obvious proper motions of the jet com- 
ponents at 2.3 and 8.4 GHz within the uncertainties. B05 also 
find no obvious superluminal motions of the components B05- 
B2 (we add "B05" before the corresponding component name 
used in their paper, this kind of nomenclature is also used for 
the components in J01 and J05), B05-B3, B05-C2, and B05-C3 
except that B05-C1 has the apparent velocity 12.1+8.0c through 
their three epochs (2003.78, 2003.83, and 2004.08) observations 
at 22 and 43 GHz. 

To show a superluminal motion of lc in 3C 66 A, a proper 
motion of at least 0.037 masyr -1 should be detected. But the 
errors of proper motions in Table |4] are all greater than 0.037 
mas yr 1 except for component b at 8.4 GHz, which means that 
the non-detection of superluminal motion may be due to the 
short time baseline of our observations, so we studied the kine- 
matics of the jet components by linking our observations with 
the observations of J01, J05, and B05. We investigated the rela- 
tions among our components and the components in B05, J01, 
and J05 by extrapolating the detected proper motions in J01 and 
J05 because of their more frequent and longer-spanned observ- 
ing epochs. Figure [4] shows the relations of these components, 
where the data points of components a and c are from our ob- 
servations at 22.2 GHz, the data points of components b and d 
are from our observations at 8.4 GHz, and other data points are 
from observations at 43 GHz in J01, J05, and B05. Although 
there is a difference in the separations from the core of the same 
component at 8.4/43 GHz and at 22.2/43 GHz probably due to 
the SSA, we show in Sect. 3.3 that such a positional difference 
is small (less than 10% of the separation of the jet component 
relative to the core), so combining the data at 8.4, 22.2, and 43 
GHz together to discuss the proper motion barely affects on our 
results. We can see from Fig. g] that J05-C4, J05-C3, J05-C2, 
J05-C1, J05-B6, J05-A1, J05-B4, and J05-A2 may be identified 
with B05-C3, a (and B05-C2), B05-C1, b, J01-B6, c, J01-B4, 
and d, respectively. We also try the linear fitting to the separa- 



tions of these components based on our identifications and give 
the results in Table 

According to Fig.@] we speculate that J01-B1 may be identi- 
fied with J05-A2. Although B05-B3 may be identified with J05- 
A2 according to the motion trend of J05-A2, B05-B2 is more 
possibly identified with J05-A2 because their position angles are 
almost the same, ~ -173°. It is also possible that J05-A1 or 
J05-C1 may be identified with B05-C1. Obviously future high- 
resolution VLBI observations are needed to confirm these iden- 
tifications. 

3.3. Frequency-dependent positional difference of the VLBI 
components 

It can be seen from Table [2] that the separations of components 
b and d from the core at 8.4 GHz are smaller than the ones at 
22.2 GHz, and the separations of component d at 2.3 GHz are 
also smaller than at 8.4 GHz; i.e., there is a difference of position 
from the core of the same jet component at different frequencies. 
There is no obvious frequency-dependent positional difference 
for the outer jet components, which may be due to their weak 
emissions and thus large positional uncertainties. According to 
Lobanov (1998), the SSA can cause the offset of the position 
from the jet apex of the ultracompact VLBI core at different fre- 
quencies due to the frequency-dependent opacity effects since 
the VLBI core is observed at the location where the optical depth 
of SSA is 1. Consequently the separations from the core of the 
optically-thin jet component at different frequencies will be dif- 
ferent. In the following, we adopt this model to estimate the sep- 
arations of 3C 66A core from the jet apex at different frequen- 
cies. 

Under the assumption that the magnetic field and electron 
density within the ultracompact VLBI core region decrease with 
r, the distance from the conical jet apex, as B oc r~ m and N oc r~" 
(Konigl ["19811 , the projected angular distance of the core from 
the jet apex on the sky plane at the observing frequency v can be 
expressed as (Lobanov |1998l l: 

l/k r l/k, 

v, v, 

r a , P mj (mas) = Ar mas (vi, v 2 )-^- k ^r- ■ v ' ' , (1) 

where k r = ((3 - 2a)m + In — 2)/(5 - 2a), and Ar mns (vi , vj) 
is the offset of the compact core position measured from the 
jet apex at frequencies v\ and V2 in the units of mas. This is 
eventually equivalent to the observed positional difference in 
the optically-thin jet component at frequencies v\ and v%, so the 
true, frequency-independent angular distance of the jet compo- 
nent from the jet apex is: 

# a (mas) =R + r a ,proj, (2) 

where R is the observed jet separation from the optically-thick 
core (Col. (4) in Table |2]i at frequency v. By using Eqs. (1) 
and (2) and adopting k r -\ in the case of the equipartition be- 
tween particle and magnetic-field energy densities, the relocated 
separations R fl of the jet component d can be calculated from 
the observational data at 2.3 and 8.4 GHz. These results found 
in Table [6] are actually the same within their errors for all six 
epochs. By substituting values of Ar mas between frequencies 2.3 
and 8.4 GHz (listed in Table [6) into Eq. (1), we can further esti- 
mate the separation of the core from the jet apex (r a ,proj) at 22.2 
GHz. These are 0.033 and 0.030 mas at two epochs 2001.20 and 
2001.48, respectively. Compared with those estimates of * a ,proj 
at 2.3 and 8.4 GHz (Tabled, we can infer from Eq. (2) that the 
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2.3 GHz 



+- tz 



2001.2 2001.4 2001.6 2001.8 2002.0 2002.2 2002.4 2002.6 

Time / yr 




2001.2 2001.4 2001.6 2001.8 2002.0 2002.2 2002.4 2002.6 
Time / yr 



Fig. 3. The linear fitting to the separations of components b, d, e, f, and g at 2.3 and 8.4 GHz, respectively. Fitting results are given 
in Table H 

Table 4. The proper motions of the jet components in 3C 66A during our six epochs. 





bs.4 


d 2 .3 


dg.4 


e2.3 


eg.4 


f 2 .3 


f 8 .4 


g2.3 


g8.4 


yu (mas yr 


') 0.038+0.029 


-0.12+0.10 


-0.048±0.040 


0.14±0.35 


0.28+0.15 


-0.10±0.57 


0.66+0.16 


0.60±1.25 


-0.06±0.61 


Papp 


1.02±0.78 


-3.21+2.68 


-1.28+1.07 


3.75±9.37 


7.50+4.02 


-2.68±15.27 


17.68±4.28 


16.07±33.48 


-1.61 + 16.34 




1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 

Time / yr 



J05-C4 

J05-C3 

J05-C2 

J05-C1 

J05-B6 

J05-A1 

J05-B4 

J05-A2 

a 

b 

d 

c 

B05-B2 
B05-B3 
B05-C1 
B05-C2 
B05-C3 
J01-B1 
J01-B4 
J01-B6 



Fig. 4. Identification of our jet components with the jet components detected in J01, J05, and B05. The open symbols represent our 
components (a, b, c, and d), filled symbols represent the components (J05-C4, J05-C3, J05-C2, J05-C1, J05-B4, J05-B6, J05-A1, 
and J05-A2) in J05, the half-filled symbols represent the components (J01-B1, J01-B4, and J01-B6) in J01, and the line-composed 
symbols represent the components (B05-B2, B05-B3, B05-C1, B05-C2, and B05-C3) in B05. See text for more details. 



Table 5. The proper motions by linearly-fitting to the separations of the jet components 
in Fig.@]by assuming the identifications of these jet components. See text for details. 
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Papp 


0.62±0.08 


2.49±0.16 


4.55+0.27 


1.85±0.27 
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-2.36±0.43 


19.28+1.07 


-2.95±0.27 
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positional differences in R are 0.29 and 0.26 mas between 2.3 
and 22.2 GHz at 2001.20 and 2001.48, and 0.054 and 0.049 mas 
between 8.4 and 22.2 GHz at 2001.20 and 2001.48, which are 
consistent with the observed differences (Ar mas ) at 2.3/22.2 GHz 
and 8.4/22.2 GHz (see Tablef2|. Similarly, with the average po- 
sitional difference of component d between 2.3 and 8.4 GHz of 
0.27+0.03 mas, we can obtain the average separation of the core 
from the jet apex (T a , P roj) at 43 GHz of 0.020 mas. So the posi- 
tional differences of the jet components from the core between 
8.4 and 43 GHz and between 22.2 and 43 GHz are only 0.08 
and 0.01 mas, which are less than 10% of the measured sepa- 
ration of component d from the core. Therefore, we can ignore 
such a frequency-dependent jet-position difference in the proper 
motion estimation discussed in Sect. 3.2. 

We must emphasise that the formula for calculating the po- 
sitional difference of the jet component is based on the assump- 
tion that the magnetic field and electron density continuously 
decrease with r; i.e., there is a smooth-flow jet structure within 
the ultracompact VLBI core region. If this is not the case, we 
still can observe the positional difference of the optically-thin jet 
component relative to the core because of the opacity effect, but 
its dependency on frequency would not follow that in Eq. (1). 

4. Discussion 

Our VLBA maps of 3C 66A show a typical core-jet structure 
with two bendings at 1 .2 and 4 mas from the core. Furthermore, 
the VLBA maps at 43 GHz in J05 show that the first bending is 
a continuous bending beginning at 1.2 mas up to 3 mas from the 
core. We also find that the flux densities of the VLBI components 
in J05, B05, and our data at 22.2 GHz decrease along the jet, but 
start to increase at the location 1.2 mas up to 2.6 mas from the 
core. Such an obvious increase in the flux densities of the VLBI 
components in J05 beginning from 1.2 mas from the core can 
be seen in Fig. [5] The coincidence between the place of the first 
bending and where the flux densities of the jet components in- 
crease suggests that the observed bending may be related to the 
increase in the flux densities of the jet. And continuous bend- 
ing from 1.2 to 3 mas may be due to the continuous decrease 
in the viewing angle to the jet components, then the increase of 
Doppler boosting. It is difficult to judge whether there is an in- 
crease in the flux density to be related to the bending at 4 mas 
since we cannot find any obvious increase in the flux density 
beginning from 4 mas. 

Combining our data at 2.3, 8.4, and 22.2 GHz at 200 1.48 
with the data at 5 GHz at 1995.65 in Taylor et al. ( 119961 1, we 
used the SSA model to fit the spectra of the core for the first 
time. The best-fit optically-thin spectral index is -0.08, and B05 
give an average spectral index -0.15 from their observations at 
22 and 43 GHz at three epochs 2003.78, 2003.83, and 2004.08. 
The difference in the spectral index may be due to the different 
core activities at different observing epochs. By examining the 
light curves from the UMRAO and the Metsahovi radio obser- 
vatory (see Terasranta et al.[T998 and Terasranta et al. 2005]), we 
find that the observing epoch 2001.48 is near the radio burst at 
2001 .20. The observations in B05 at three epochs were also close 
to a radio burst at 2004.0, but the amplitude of this radio burst is 
lower than the one at 2001.20. It may be possible that there are 
more high-energy electrons in the core region at 2001.48 than at 
three epochs in B05, resulting in a flatter spectrum at 2001.48 
than those at 2003.78, 2003.83, and 2004.08. 

We fit the separations of the jet components of 3C 66A at 
2.3 and 8.4 GHz for six epochs, and did not detect any proper 
motions of these jet components during our observations, which 



may be due to our short coverage time. We then studied the re- 
lations between our jet components and the components in J01, 
J05, and B05, and give the possible identifications of these jet 
components in Fig. [4] We can see from Fig.|4]that our jet com- 
ponents are not stationary if our identifications are right. The 
components J05-A1 and J05-A2 in J05 show the apparently in- 
ward motion to the core, it may be due to the viewing angle 
to the jet component becoming small when the jet component 
moves back toward the line of sight. This is also evidence that 
the increase in the flux densities of the jet components from 1.2 
to 2.6 mas may be due to decreasing viewing angles to these jet 
components. 

The average positional difference in component d from the 
core at 2.3 and 8.4 GHz is 0.27 mas. Lobanov d 19981 ) give the 
positional differences of the optically-thin jet components from 
the core at 2.3 and 8.4 GHz of 3C 395 (z=0.635), 4C 39.25 
(z=0.699), and 1038+528A (z=0.678) about 0.60, 0.50, and 0.50 
mas, respectively, which are bigger than the one, 0.27 mas, in 
3C 66A. The value of the observed positional difference of the 
jet component from the core is determined by the magnetic field 
and electron density in the core region, redshift, and viewing 
angle to the jet. The observed positional difference of the jet 
component will be bigger with the stronger magnetic field and 
higher electron density in the core region, the lower redshift, and 
the larger viewing angle. The redshift of 3C 66A is lower than 
3C 395, 4C 39.25, and 1038+528A, but the magnetic field in the 
core region and the viewing angle are smaller than the ones in 
3C 395, 4C 39.25, and 1038+528A (Lobanov [T99"8l ), which may 
explain why the observed positional difference in the component 
d in 3C 66A is smaller than the ones in 3C 395, 4C 39.25, and 
1038+528A. 

5. Conclusions 

We present results of VLBA observations at 2.3, 8.4, and 22.2 
GHz during six epochs (the observations at 22.2 GHz were only 
performed in the first two epochs). We find from our VLBA 
maps that there are two bendings at 1 .2 and 4 mas from the core. 
We identify the brightest k component with the core from our 
simultaneous spectral data that show an inverse. Combining the 
existing VLBI data at 5 GHz, we use the SSA model to fit the 
core spectrum at epoch 2001.48 for the first time, and obtain the 
best-fit optically-thin spectral index -0.08, the magnetic filed in 
the core region 1.1 mG, and the maximum flux density 0.66 Jy 
at the turnover frequency of 6.52 GHz. We did not detect any 
obvious proper motions in jet components during our observing 
epochs, which may be due to the short observing time. We give 
a possible identification of our jet components with the jet com- 
ponents in J01, J05, and B05 by analysing the proper motions of 
these jet components. The positional differences in component d 
from the core at different frequencies are consistently seen at all 
the six epochs. 
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-173.5+0.1 


0.56+0.01 


1.0 





9.39 x 10 9 




e 


0.043+0.004 


6.39+0.32 


170.5±0.6 


3.72+0.14 


1.0 





7.76 x 10 7 




f 


0.037+0.004 


13.21+0.66 


165.3±0.5 


5.02+0.29 


1.0 





3.66 x 10 7 




g 


0.026+0.0026 


23.0+1.15 


165.6+1.5 


10.52+1.4 


1.0 





5.86 x 10 6 


Epoch 2002.45, X z vXi = 1-06, tf M = L05 


2.3 


k 


0.429+0.043 








0.30+0.05 


1.0 





1.59 x 10 12 




d 


0.209+0.021 


1.98+0.03 


-174.9±0.2 


0.81+0.04 


1.0 





1.06 x 10 11 




e 


0.054+0.005 


6.47+0.32 


168.6+0.3 


2.90+0.08 


1.0 





2.14 x 10 9 




f 


0.056+0.006 


13.14+0.66 


165.3+0.3 


5.32+0.13 


1.0 





6.59 x 10 8 




g 


0.031+0.0031 


22.1 + 1.10 


165.6+0.6 


10.12+0.45 


1.0 





1.01 x 10 8 


8.4 


k 


0.447+0.045 








0.20+0.01 


1.0 





2.79 x 10 11 




b 


0.136+0.014 


0.79+0.02 


-165.9±0.2 


0.37+0.01 


1.0 





2.48 x 10 10 




d 


0.123+0.012 


2.22+0.02 


-174.1±0.2 


0.75+0.01 


1.0 





5.46 x 10 9 




e 


0.046+0.005 


6.62+0.33 


168.4±0.5 


4.21+0.14 


1.0 





6.48 x 10 7 




f 


0.045+0.004 


13.3+0.66 


164.9±0.4 


5.27+0.21 


1.0 





4.04 x 10 7 




g 


0.021+0.0021 


22.96+1.15 


168.7+0.8 


7.07+1.78 


1.0 





1.05 x 10 7 



Notes:(l) Observing frequency; (2) component name; (3) flux density of component;(4) separation from the core (component 
k); (5) P. A. (from noth to east); (6) major axis of the elliptical Gaussian component; (7) ratio of the minor to major axes; (8) position 
angle of the major axis of the elliptical Gaussian component from north to east; (9) brightness temperature according to the formula 
given in Shen et al. ( 11997b . 



